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It was recently shown that the vacuum in the background of a strong enough magnetic field may 
O | become an electromagnetic superconductor due to interplay between strong and electromagnetic 

O | forces. The superconducting ground state of the QCDxQED sector of the vacuum is associated 

with magnetic-field-assisted emergence of quark-antiquark condensates which carry quantum 
i— i numbers of charged p mesons (i.e., of electrically charged vector particles made of lightest, u and 

d, quarks and antiquarks). Here we demonstrate that this exotic electromagnetic superconduc- 
tivity of vacuum is also accompanied by even more exotic superfluidity of the neutral p mesons. 
The superfluid component - despite being electrically neutral - turns out to be sensitive to an 
external electric field as the superfluid may ballistically be accelerated by a test background elec- 

CN 

trie field along the magnetic-field axis. In the ground state both superconducting and superfluid 
components are inhomogeneous periodic functions of the transversal (with respect to the axis of 
CNI the magnetic field) spatial coordinates. The superconducting part of the ground state resembles 

an Abrikosov ground state in a type-II superconductor: the superconducting condensate orga- 
nizes itself in periodic structure which possesses the symmetry of an equilateral triangular lattice. 
Each elementary lattice cell contains a stringlike topological defect (superconductor vortex) in 
the charged p condensates as well as three superfluid vortices and three superfluid antivortices 
made of the neutral p condensate. The superposition of the superconductor and superfluid vortex 
lattices has a complicated "kaleidoscopic" pattern. 
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1. Introduction 

Properties of Quantum Chromodynamics in an extremely strong external magnetic field of 
hadronic scale have attracted strong interest of the scientific community. For example, hot quark- 
gluon plasma created in heavy-ion collisions may exhibit the chiral magnetic effect [1]: chirally- 
imbalanced matter may generate electric current along the axis of the magnetic field [2]. The 
attractive feature of the chiral magnetic effect is that it may be tested experimentally: the magnetic 
field of high strength appears naturally in noncentral collisions of heavy ions while the chiral im- 
balance is generated by topological transitions which are inherent to QCD [1, 3]. Moreover, the 
phase diagram of QCD matter [4, 5] is sensitive to the presence of strong magnetic fields which 
opens a possibility to check these effects experimentally in the astrophysical setup (the magnetic 
field affects the macroscopic properties of the neutral stars, such as mass, adiabatic index, moment 
of inertia, and cooling curves [5]). 

The magnetic field has an impact not only on the properties of the (quark) matter, but it may 
also affect the quantum vacuum. For example, strong magnetic field enhances the chiral symmetry 
breaking due to magnetic catalysis phenomenon [6, 7, 8]. Consequently, the background magnetic 
field changes the finite-temperature phase structure of the QCD vacuum as it shifts the critical 
temperatures of the confinement-deconfinement and chiral transitions 1 [9, 10, 11, 12, 13]. 

Another interesting phenomenon is that the sufficiently strong magnetic field may induce a 
spontaneous electromagnetic superconductivity of the vacuum [14, 15]. The superconductivity 
of, basically, empty space, is supported by emergence of specific quark-antiquark condensates 
with the quantum numbers of p mesons (electrically charged vector particles made of lightest, u 
and d, quarks and antiquarks). The spontaneous creation of both charged and neutral p-meson 
condensates should take place if the strength of the magnetic field exceeds the critical value 

fi f ~10 16 Tesla or efi c ~ 0.6GeV 2 . (1.1) 

It is very promising that the magnetic fields of (the twice of) the required strength (1.1) may be 
created in the experiments on heavy-ion collisions at the Large Hadron Collider [16]. 

Below we discuss the complicated structure of the superconducting ground state in the scope 
of the effective ("bosonic") quantum electrodynamics for the charged and neutral p mesons [17] 
following Refs. [14, 18]. Similar results can also be obtained in a completely different ("fermionic') 
approach which uses the Nambu-Jona-Lasinio model [15]. This somewhat counterintuitive effect 
("nothing" becomes a superconductor) can also be found in holographic approaches [19, 20] and 
in numerical simulations of quenched QCD on the lattice [21]. 

Interestingly, the ground state contains also a neutral superfluid component [22]. This super- 
fluid turns out to have a rather exotic property as it can be accelerated by an external electric field. 

Another exotic property of this ground state was pointed out recently by Smolyaninov: the 
superconducting state of the vacuum behaves similarly to a (hyperbolic) metamaterial which may 
act as diffractionless "perfect lenses" [23]. Since large magnetic fields may have existed in the 
early moments of our Universe [24], a possible metamaterial behavior of the early vacuum may 
have certain imprints in the large scale structure of the present-day Universe. 

'One should notice, however, that up to now there is no general consensus in the literature even on a qualitative 
behavior of these transitions in the background of the strong magnetic field. 
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2. Electrodynamics of charged and neutral p mesons 

The p meson is a vector (spin-1) particle made of a light (u or d) quark and a light anti- 
quark. There are positively, p + = ud and negatively, p~ = du charged p mesons and a neutral 
dd)/V2 meson with approximately the same masses m p rj 775 MeV. A self-consistent 
quantum electrodynamics for these charged and neutral vector mesons is described by Djukanovic- 
Schindler-Gegelia-Scherer (DSGS) Lagrangian [17]: 

* = -\ w v - \ pUp^+K P y - \ P <V 0) * V + f p£V 0)m + ^ f^ p (? , (2.D 

which generalizes the celebrated vector meson dominance model [25] by including the electromag- 
netic (Maxwellian) t/(l) e . m . group: 

.(.).., {>f^f:^f/ M (2.2) 

The Lagrangian (2. 1) describes the interaction of the electromagnetic field with the neg- 
atively charged, p M = p~ = (p^, 1 ' — ipjj/)fy/2, positively charged, p+ = p^, and neutral, pf , 
vector mesons, as well as the self-interactions of the mesons. The tensor quantities in Eq. (2.1) are 
various strength tensors, 

Fptv = <Vv - d v A^ , f$ = d M pf - d v p£ 0) , (2.3) 

p£v = /£v - igsipjiPv - PuPv) > Pmv = d mPv - -DvPm , (2.4) 

and the covariant derivative is = + ig s p^ — ieA^. The p7T7T vertex coupling g s is fixed to its 
phenomenological value g s = g pjtK « 5.88. 



3. Magnetic-field-induced condensation of charged p mesons 

There is a very simple naive argument which supports the condensation of the p meson ex- 
citations [14]. Consider a free charged relativistic vector (spin-1) particle. Assume that it has the 
mass m and the gyromagnetic ratio g = 2 (which is strongly supported by various results [14, 18]) . 
Then in a background of an external magnetic field B ext = (0,0,B ext ) this particle has the following 
energy spectrum: 

e« A {Pz) = P\ + ( 2 « - 2s z + 1 > fi ext + m 2 . (3.1) 

where s z = — 1,0, +1 is the spin projection on the field's axis z = %3, n ^ is a nonnegative integer 
number (which, together with the spin projection s z , labels the Landau levels), and p z is the particle 
momentum along the magnetic field's axis £. The ground state (given by the quantum numbers 
n = 0, s z = +1, p z = 0) has the following (squared) energy: £ 2 (S ext ) = m 2 — eB ext . The ground 
energy of the charged p -meson in the external magnetic field becomes purely imaginary quantity 
if the magnetic field exceeds the critical value B c = m 2 je, which, in the case of the p mesons, 
corresponds to the critical field (1.1). Exactly the same result can be derived from the DSGS 
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model (2.1) in the background of the magnetic field 2 [14]. The condensed p -meson state has the 
following wavef unction: 

Pi(xi,x 2 ) = -ip 2 (x 1 ,x 2 ) = p(xi,x 2 ) , po(xi,x 2 ) =Pi{x\,x 2 ) = 0, (3.2) 

which can be expressed via a single complex scalar function p = p(x\,x 2 ). Since the global ro- 
tation around the magnetic field axis as well as local electromagnetic gauge transformations (2.2) 
rotate the scalar field p by a phase, the p-meson condensate (3.2) breaks the gauge symmetry(2.2) 
spontaneously by "locking" the subgroup of the Lorentz rotations with the gauge transformations 
("Lorentz gauge locking" [14]). 

Notice that the similar condensation effects exist in other field theories: it is suggested to 
happen in the standard model of electroweak interactions in a background [with a much stronger 
than (1.1)] magnetic field [26, 27] . Another example is given by a pure Yang-Mills theory, in which 
the gluons should form a condensate in a c/iromomagnetic field background [28]. 

4. Superconducting ground state: an explicit solution 

The effective bosonic model (2. 1 ) in QCD plays the same role as the Ginzburg-Landau model [29] 
plays in the superconductivity: both bosonic models describe dynamics of the superconducting 
carriers (the p-meson excitations in QCD and the Cooper pairs in ordinary superconductivity). 
Following the standard approach to superconductivity, we solve the p-meson DSGS model (2.1) 
classically in order to find basic features of the electromagnetic superconductivity of the QCD vac- 
uum. In particular, we utilize the Abrikosov solution in a type-II superconductor [29], and we 
choose a general solution of classical equations of motion of the model (2.1) in a form of a sum 
over lowest Landau levels: 

P(z) = LCA/v,-^,-M , h n (v,z,z) =exp{-^(|z| 2 + z 2 ) -nv 2 n 2 + 2nvnz\, (4.1) 
neZ ^ L B L B J V 2 J 

where Lb = y / 2n/(eB) is the magnetic length and z = x\ + ix 2 is a complex coordinate in the 
transversal (with respect to the direction of the magnetic field axis, z = X3). Here v and C n are 
certain parameters which are chosen to minimize the corresponding energy functional. A symmetry 
pattern of the C n coefficients, C,, + n = C n , determines the type of the Abrikosov lattice. The solution 
is generally periodic in the {x\,x 2 ) plane and it corresponds to a certain lattice symmetry: N = 1 
gives us a square lattice, N = 2 corresponds to an equilateral triangular lattice etc. 

It turns out that in QCD - contrary to the Ginzburg-Landau model - the coefficients C n are 
dependent on the value of the magnetic field, C„ = C n (B). In Ref. [18] it was shown that close to 
the phase transition, B ^ B c , the solution is given by the equilateral triangular lattice (N = 2 with 
Ci = iC 2 and v = \/3/ y/2). Each elementary lattice cell contains a stringlike topological defect 
(superconductor vortex) in the charged p condensates as well as three superfluid vortices and three 
superfluid anti vortices made of the neutral p(°) condensate. The superposition of the superconduc- 
tor and superfluid vortex lattices has a complicated "kaleidoscopic" pattern. The behavior of the 
charged and neutral condensates in the ground state at B = 1.01B C are shown in Fig. l(left) and 
1 (right), respectively. Further details can be found in Ref. [18]. 

2 In the DSGS Lagrangian (2.1) the anomalously large gyromagnetic ratio (g = 2) of the charged p ± mesons is 
described by a nonminimal coupling of the p mesons to the electromagnetic field given by the last term in Eq. (2.1). 
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p(xi,x 2 ) p (0 >(x u x 2 ) 



Figure 1: The charged (left) and neutral (right) p-meson condensates in the ground state at B = 1 .0\B C . The 
3D plots show the behavior of the absolute values (in MeV) of the condensates while the projection on the 
(xi,x 2 ) plane is a density plot of the phases of the corresponding condensates. The white lines are gauge- 
dependent singularities (Dirac sheets) attached to the superconductor and superfluid vortices, respectively. 



The electric current density, 

J, = ie[p^p v ,-p v p^ + d v (plp,-pj l p v )]--d v fi°l (4.2) 

satisfies a London-like equation in the longitudinal (parallel to the magnetic field) direction [14]: 

dJ^(x) dJn(x) , , . . , 9 / ntn , 

- 1 p- L -^- L = -K(x h x 2 )E 3 (x), K(x 1 ,x 2 )=4e 2 ( 2 p 2 )(x h x 2 ), (4.3) 
oxq ox?, \ —dj_ + m.Q / 

where = df + d 2 is the transversal Laplacian. Relations (4.3) constitute the London-type equa- 
tion which is typical for superconducting systems in which the electric DC current flows without 
resistance. The quantity K* = k(x\ ,x 2 ) is a nonlocal function of the superconducting condensate. In 
the transversal directions the ground state behaves as an insulator. 

The condensation of the charged (superconducting) field p^ induces condensation of a neutral, 
superfluid-like field p^ in the ground state of the vacuum [14]: 

p w {x u x 2 )=pf )) (x u x 2 ) + ip^ ) (xi,x 2 ) = 2ig s -(—^- — 2lP| 2 )(*b*2), (4.4) 

The longitudinal components of the neutral condensate are zero, p^ = p^ = 0. 

It is interesting to note that the electrically neutral current of p(°) mesons, defined via the 
relation Jfc' = —^d y fl° l }, satisfies a London-like equation as well: 

dJ 3 °\x) djf(x)_ .^0)^^ K i*) {xuX2 )=^( J p| 2 W;c 2 ).(4.5) 



oxo a xt, \—dj_+niQ 

Literally, equations (4.3) and (4.5) indicate that if we impose a weak (test) external electric 
field parallel to the axis of the strong magnetic field, then the charged and neutral currents, respec- 
tively, would start to flow frictionlessly along the magnetic field axis. Once the electric field is 
removed the currents continue to flow because of the absence of the dissipation for the supercon- 
ducting and superfluid condensates. The sensitivity of the neutral condensate to the electric current 
is an unexpected feature of this new vacuum ground state. 
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